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We use (111 ± 1) million T{?,S) and (89 ± 1) million T{2S) events recorded by the BAa4R de- 
tector at the PEP-II B-factory at SLAC to perform a study of radiative transitions between bot- 
tomonium states using photons that have been converted to e+e" pairs by the detector material. 
We observe T(3S') 7Xf)0,2(lf') decay, make precise measurements of the branching fractions for 
Xi)i,2(lf, 2P) — > 7T(15') and X6i,2(2P) — >• ^T{2S) decays, and search for radiative decay to the 
rjb{lS) and rn,{2S) states. 

PACS numbers: 13.20.Gd, 14.40.Pq 



I. INTRODUCTION 



Bottomonium spectroscopy and radiative transitions 
between bb states can be well-described by effective po- 
tential models [1]. To leading order, radiative decays 
are expected to be dominantly electric (El) or magnetic 
(Ml) dipole transitions. In the non-relativistic limit, the- 
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oretical predictions for these decays are straightforward 
and well-understood. However, there arc a few notable 
cases where the non-relativistic decay rates are small or 
zero, e.g. in "hindered" Ml transitions between S-wave 
bottomonium such as T(nS) — !> "frii,{n' S) [n > n'), and 
as a consequence of small initial- and final-state wave- 
function overlap in the case of T{3S) — > 7X;,j(lP) decays 
[2] ; higher-order relativistic and model-dependent correc- 
tions then play a substantial role. Measurements of these 
and other El transition rates can lead to a better under- 
standing of the relativistic contributions to, and model 
dependencies of, interquark potentials. Furthermore, be- 
cause radiative transitions have a distinct photon energy 
signature associated with the mass difference between the 
relevant bb states, they are useful in spectroscopic studies 
for mass measurements, and in the search for and iden- 
tification of undiscovered resonances. 
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Radiative transitions within the bottomonium system 
have been studied previously in several experiments, sueh 
as Crystal Ball [3, 4], ARGUS with converted photons [5], 
and iterations of CUSB [6-10] and CLEO [11-15] (includ- 
ing an analysis of photon pair conversions in a lead ra- 
diator inserted specifically for that purpose [16]). These 
analyses have focused mainly on Xbj(nP)-related mea- 
surements, such as the determination of the masses and 
the El transition rates to and from T{mS) states. More 
recently, the BaBAR experiment finished its operation by 
collecting large samples of data at the Y{3S) and T{2S) 
center-of-mass (CM) energies. These data are useful for 
studies of bottomonium spectroscopy and decay and have 
already led to the discovery of the long-sought 7?&(1<S') 
bottomonium ground state [17, 18], an observation later 
confirmed by CLEO [19]. 

In this paper, we present a study of radiative transi- 
tions in the bottomonium system using the inclusive con- 
verted photon energy spectrum from T(35) and T{2S) 
decays. The rate of photon conversion and the recon- 
struction of the resulting e+e~ pairs has a much lower 
detection efficiency than that for photons in the BaBAR 
electromagnetic calorimeter, a disadvantage offset by 
a substantial improvement in the photon energy reso- 
lution. This improvement in resolution is well-suited 
for performing precise transition energy (hence, parti- 
cle mass, and potentially width) measurements, and to 
disentangle overlapping photon energy lines in the in- 
clusive photon energy spectrum. This analysis has dif- 
ferent techniques, data selection, and systematic uncer- 
tainties than the previous studies [17 19], and is rela- 
tively free from complications due to overlapping tran- 
sition peaks, and calorimeter energy scale and mea- 
surement uncertainties. We report measurements of 
Xbj{2P) ^ -rT{2S), Xbj{iP,2P) ^ 7^(1^), observation 
of r{3S) 7X60,2(lf), and searches for the r]bilS,2S) 
states. 

In Sec. II we describe the BaBAR detector and the 
data samples used in this analysis. Section III describes 
the photon conversion reconstruction procedure and the 
event selection criteria. Each of the following sections 
(Sec. IV - VII) individually describes the analysis of a 
particular region of interest in the inclusive photon en- 
ergy spectrum. Section VIII summarizes the results ob- 
tained. Appendix provides specific details of some sys- 
tematic uncertainties related to this analysis. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

The BaBar detector is described in detail elsewhere 
[20] ; a brief summary is provided here. Moving outwards 
from the collision axis, the detector consists of a double- 
sided five-layer silicon vertex tracker (SVT) for measur- 
ing decay vertices close to the interaction point, a 40- 
layer drift chamber (DCH) for charged-particle tracking 
and momentum measurement, a ring-imaging Cherenkov 



detector for particle identification, and a CsI(Tl) crystal 
electromagnetic calorimeter (EMC) for measuring the en- 
ergy deposited by electrons and photons. These detector 
subsystems are contained within a large solenoidal mag- 
net which generates a 1.5-T field. The steel magnetic 
flux return is instriimented with a muon detection sys- 
tem consisting of resistive plate chambers and limited 
streamer tubes [21]. 

The inner tracking region also contains non- 
instrumented support structure elements. Interior to the 
SVT, the interaction region is surrounded by a water- 
cooled, gold-coated beryllium beam pipe. The SVT 
support structure consists primarily of carbon-fiber and 
Kevlar®. The SVT, beam pipe and vacuum cham- 
ber, and the near-interaction-point magnetic elements 
are mounted inside a cylindrical, carbon-fiber support 
tube. The inner wall of the DCH is a cylindrical tube of 
beryllium coated with anti-corrosion paint. A photon at 
normal incidence traverses approximately 0.01 radiation 
lengths (^o) of material before reaching the SVT, and 
an additional Q.OSXq before the DCH. Due to the asym- 
metric energy of the incoming e+e~ beams, the photons 
in this analysis tend to be boosted in the direction of e~ 
beam, increasing the typical number of radiation lengths 
up to 0.02Xo and O.OSXq to reach the previously noted 
detector subsystems. While this extra material is usually 
considered detrimental to detector performance, it is es- 
sential for 7 e+e~ conversions in the present analysis. 

The BaBAR detector collected data samples of (121 ± 1) 
million r{3S) and (98 ±1) million T{2S) decays [22] pro- 
duced by the PEP-II asymmetric energy e+e~ collider. 
This corresponds to an integrated luminosity of 27.9±0.2 
fb-i (13.6 ± 0.1 fb-i) taken at the r{3S) {r{2S)) res- 
onance. Approximately 10% of these data (referred to 
here as the "test sample") were used for feasibility stud- 
ies and event selection optimization; they are excluded in 
the final analysis. The results presented in this analysis 
are based on data samples of (111 ± 1) million T{3S) and 
(89 ± 1) million r{2S) decays. An additional 2.60 ± 0.02 
(1.42 ±0.01) fb^^ of data were taken at a CM energy ap- 
proximately 30 MeV below the nominal T{3S) {T{2S)) 
resonance energy, to be used for efficiency-related studies. 

Large Monte Carlo (MC) datasets simulating the sig- 
nal and expected background decay modes are used for 
the determination of efficiencies and the parameteriza- 
tion of lineshapes for signal extraction. The particle pro- 
duction and decays are simulated using a combination of 
EVTGEN [23] and JETSET [24]. The radiative decays 
involving Xbj{nP) states are assumed to be dominantly 
El radiative transitions, and the MC events are gener- 
ated with theoretically predicted helicity amplitudes [25]. 
The interactions of the decay products traversing the de- 
tector are modeled by Geant4 [26]. 
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III. EVENT RECONSTRUCTION AND 
SELECTION 

Photon conversions are reconstructed with a dedicated 
fitting algorithm that pairs oppositely charged particle 
tracks to form secondary vertices away from the interac- 
tion point. The algorithm minimizes a value (x/it) 
based on the difference between the measured helical 
track parameters and those expected for the hypothe- 
sis that the secondary vertex had originated from two 
nearly parallel tracks emitted from a 7 — >■ e+e^ conver- 
sion. The Xfit value includes a term to account for an ob- 
served finite opening angle between the converted tracks. 
Requiring xju < 34 is found to be the optimal value to 
select a high-purity converted photon sample. The re- 
constructed converted photons are also required to have 
an e+e^ invariant mass of mg+g- < 30MeV/c^ (though 
in practice, mg+e- is typically less than lOMeV/c^). To 
remove internal conversions and Dalitz decays, and to 
improve signal purity, the conversion vertex radius (pj) 
is required to satisfy 1.7 < pj < 27 cm. This restricts 
the photon conversions to the beampipe, SVT, support 
tube, and inner wall of the DCH, as seen in the plot of 
conversion vertex position for a portion of the "test sam- 
ple" in Fig. 1. The efficiency for photon conversion and 
reconstruction versus energy in the CM frame (E*), as 
determined from a generic Y(SS) MC sample, is shown 
in Fig. 2. 

Figure 3 shows the inclusive distributions of the re- 
sulting reconstructed converted photon energy. The data 
are divided into four energy ranges, as indicated by the 
shaded regions in Fig. 3. These ranges and the cor- 
responding bottomonium transitions of interest are, in 
T{3S) data: 

• 180 <E*< 300 MeV: Xb.j{2P) -/T{2S) 

• 300 < E* < 600 MeV: T{3S) lXbj{lP) and 
T(3^)^77/h(25) 

• 600 < E* < 1100 MeV: X6j(2-P) ^ lT{lS) and 
T(35) ^ iVb{lS) 

and in T{2S) data: 

• 300 < E* < 800 MeV: XbJ^^P) ^ 72"(15') and 
T(25)^ 777fc(15). 

Figure 4 summarizes these energy ranges and the radia- 
tive transitions of interest in a pictoral form. Peaks re- 
lated to some of these transitions are already clearly vis- 
ible in Fig. 3, where the photon energy in the CM frame 
of the initial particle for the radiative transition from 
an initial (i) to final (/) state is given in terms of their 
respective masses by 

9 9 
m, — rri f „ 

E,{i^f)^^ U\ (1) 

2mi 

Because we analyse the photon energy in the CM frame of 
the initial T{mS) system {E*), the photon spectra from 




■^So -20 -10 10 20 30 ^ 

X (cm) 

FIG. 1: (Color online) End view of the the BABAR inner detec- 
tor along the beam axis as seen by converted photons. Points 
indicate the number of converted photon vertices per cross- 
sectional measured in a subset of the "test sample" 
data. From the center outwards, features of note include the 
beam pipe, the SVT {e.g. hexagonal inner layers) and its sup- 
port structure rods, the support tube, and the inner wall of 
the DCH. 



subsequent boosted decays {e.g. Xbj{nP) — > 7^(15)) are 
affected by Doppler broadening due to the motion of the 
parent state in the CM frame. 

To best enhance the number of signal {S) to back- 
ground {B) events, the event selection criteria are cho- 
sen by optimizing the figure of merit J- = --y==. 
This is done separately for each energy region. The 
180 < E* < 300 MeV energy region in r{3S) uses the 
same criteria as determined for the similarly low energy 
300 < E* < 600 MeV range. We determine S from MC 
samples of T{mS) '-fr]b{lS) weighted to match the 
measured branching fractions [27], and T{3S) — >■ jr]b{2S) 
assuming the same branching fraction as for the decay to 
jr]h{lS). Because the generic decay processes of T{3S) 
are not well-known (ie: a large percentage of the exclusive 
branching fractions have not been measured), the "test 
sample" data are used to estimate B. The optimization 
is performed by varying the selection criteria for the to- 
tal number of tracks in the event {nTRK), the absolute 
value of the cosine of the angle in the CM frame between 
the photon momentum and the thrust axis (|cos0t|) 
[28] , and a 7r° veto excluding converted photons produc- 
ing an invariant mass {ruj-y) consistent with m^o when 
paired with any other photon (converted or calorimeter- 
detected) above a minimum energy {E^2) in the event. 
A requirement on the ratio of the second and zeroth Fox- 
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FIG. 2: Efficiency for tlie conversion and reconstruction of 

a plioton versus photon energy as derived from a sample of 
generic T{3S) MC events before the optimal selection criteria 
have been applied. 




FIG. 4: Pictoral representation of energy levels in the bot- 
tomonium system and the radiative transitions studied here. 
The Roman numeral labels indicate the corresponding Section 
in which the given transition is discussed. 




TABLE I: Acceptance criteria for converted photon events. 



Variable 


e; 


Range (MeV) 




r(35) 


r{3S) 


r{2S) 




[180,600] [600,1100] [300,800] 


nTRK 


> 8 


> 8 


> 8 


\ cos 9t \ 


< 0.85 


< 0.75 


< 0.85 


\m^^ - m^o| (MeV/c^) 


> 10 


> 20 


> 20 


E^2 (MeV) 


> 90 


> 75 


> 70 


R2 


< 0.98 


< 0.98 


< 0.98 




Ej(GeV) 



FIG. 3: Raw inclusive converted photon energy spectrum 
from (a) Y{3S) and (b) T(2S) decays. The shaded areas in- 
dicate different regions of interest considered in detail in this 
analysis. The Roman numeral labels indicate the correspond- 
ing Section in which each energy region is discussed. 



Wolfram moments [29] of each event, i?2, is also applied. 
The reason for using these particular variables (indicated 
in parentheses) is to preferentially select bottomonium 
decays to hadronic final states (nTRK) and to remove 
photons from continuum background events (| cos 9t \ and 
R2) and 7r° decays (m^o veto). Table I summarizes the 
values for the optimized selection criteria. 

The efficiency for reconstruction and selection of signal 
events (e) is determined from MC simulation. A dedi- 
cated e'^e~ — > ii'^fi^j sample is used to study our de- 
tector model and converted photon efficiency (discussed 
in Appendix A), and the correspondence between simu- 
lation and data is found to be in very good agreement. 
Once the optimal selection criteria have been applied, 
e < 1.5% for conversions compared to ~ 40% for photons 
in the EMC for the energy range of interest in this analy- 
sis. Conversely, a large improvement is gained in photon 
energy resolution, e.g. from ~ 25 MeV in the calorime- 
ter to 4 MeV or better with converted photons. Figure 
3 demonstrates both of these features. The sharply- 
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peaking structures correspond to bottomonium transi- 
tions, and are narrow and well-resolved in this analy- 
sis. Unlike in the photon energy spectrum expected from 
the EMC [17, 18], the distribution for converted pho- 
tons drops with energy. The efficiency decreases (also 
seen in Fig. 2) due to the inability to fully reconstruct 
the conversion pair as at least one of the individual track 
momenta approaches the limit of detector sensitivity. We 
are unable to contribute useful new information on tran- 
sitions expected below E* =180 (300) MeV for the T{3S) 
(T(25)) analysis, which is why those energy ranges are 
not considered here. 

The number of signal events for a given bottomonium 
transition is extracted from the data by performing a 
fit to the E* distribution in 1 MeV bins. The functional 
form and parameterization for each photon signal is de- 
termined from MC samples, as described below. In gen- 
eral, the lineshape is related to the Crystal Ball function 
[30], i.e. a Gaussian function with a power-law tail. This 
functional form is used to account for bremsstrahlung 
losses of the e+e~ pair. Comparisons between simula- 
tion and data made on e+e~ — > e+e~ events used for the 
standard luminosity measurement in BaBAR demonstrate 
that the bremsstrahlung tails of these distributions are 
found to be well-described. The underlying smooth in- 
clusive photon backgroimd is described by a fourth-order 
polynominal multiplied by an exponential function. This 
functional form adequately describes the background in 
each separate energy range. 



IV. r(3S) : 180 < -E; < 300 MeV 

The main purpose of the fit to the 180 < E* < 
300 MeV region of the T{3S) photon energy spectrum, 
shown in detail in Fig. 5, is to measure the Xb.ii'^P) 
7T(2S') transitions. The only previous measurements of 
these transitions were made by CUSB [10] and CLEO 
[12] nearly two decades ago. Those analyses examined 
the low-energy photon spectrum from exclusive T{3S) — > 
77T(2S')(£+£~) decays to derive the branching fractions 
for B(x6i,2(2P) 7T(2S')), and in the case of the CUSB 
result, to obtain evidence for Xba{2P) jT{2S). We 
present the first fit to E* to measure the photon from 
X6j(2P) — >■ 7T(2S') directly. Though this analysis is po- 
tentially sensitive to all six T{lDj) — > 7Xtj(l-P) decays, 
we treat these decays as a small systematic effect to the 
Xbj(2P) -T- 7T(2S') measurement. 

The X6j(2P) transition lineshapes are parameterized 
by a Gaussian with power law tails on both the high and 
low side. This is best understood as a "double-sided" 
Crystal Ball function with different transition points and 
exponents for the high and low tails, but with a common 
Gaussian mean and standard deviation in the central re- 
gion. The effects of Doppler broadening, due to the mo- 
tion of the X6j(2P) in the CM frame, are small (~ 2 MeV 
width) for these transitions. The T(ir'j)-related line- 
shapes are individually parameterized in terms of a single 



Crystal Ball function. Parameterization of these transi- 
tions presents a complication because only the mass of 
the J = 2 state has been measured reliably [31, 32], the 
value mr(iD) = (10163.7 ± 1.4) MeV/c^ being obtained 
when the experimental results are averaged. Marginal 
evidence for the J = 1 and 3 states was also seen 
at ~ 10152 MeV/c^ and ~ 10173 MeV/c^, respectively 
[31, 32]. These values are consistent with several the- 
oretical predictions [33]. given a shift to bring the the- 
oretical value for mx(\D^) into agreement with experi- 
ment. We therefore assume the mn^iDi 3) mass values 
stated above to compute the expected energy for transi- 
tions from those states. The event yields for these transi- 
tions are fixed to the branching fractions expected when 
S(r(3S') 7Xfcj(2P)) [27] is combined with the pre- 
dictions for B{xbj{2P) ^ 77X6,7(1P)) via r{lDj) [34]. 
The efficiencies for the T{1D) transition signals range 
from approximately 0.17 to 0.30%, monotonically rising 
with E*. 

Figure 5 shows the measured photon spectrum and re- 
sults of the fit, before and after subtraction of the inclu- 
sive background. In this fit, the parameters describing 
the background and any systematic offset in the; E* scale 
are free parameters, together with the signal yields for 
X6,/(2P) 7T(25') decays. Table II summarizes the 
fit results. Considering both statistical and systematic 
uncertainties, we find significant Xbi,2(2P) — >• jT{2S) 
signals (> 12a and > 8a, respectively, where a repre- 
sents standard deviation), but do not find evidence for 
Xbo(2P) — > 7^(25') decay. The overall energy offset, de- 
termined predominantly by the position of the X6i,2(2P) 
transition peaks compared to the nominal [27] values, is 
found to be inconsequential (— 0.3 ± 0.2 MeV). 

The systematic uncertainties on these measurements 

(with their approximate sizes given in parentheses below 
and throughout) include the uncertainty in the fit param- 
eters fixed from MC, uncertainty in the converted pho- 
ton efficiency, assumptions related to the T{lDj) con- 
tributions, uncertainty on masses used to calculate the 
expected E* values, the T{mS) counting uncertainty, ef- 
fects of the fit mechanics, and the effect of the choice 
for the background shape. For each fit component, all 
of the parameters fixed to MC-dctermined values are 
varied individually by ±lcr of the statistical uncertainty 
from the MC determination, and the fit repeated. The 
maximal variation of the fit result for each component 
is taken as the systematic uncertainty, and summed in 
quadrature (~ 4%). The systematic uncertainty on the 
converted photon efficiency (4.7%) is estimated using an 
off-peak control sample and varied selection criteria, as 
described for all energy regions in Appendix A. The fits 
are repeated with the T{lDj) masses individually varied 
by their approximate experimental uncertainties (±1.8, 
±1.4, and ±1.5 McV/c^ for J = 1,2 and 3, respectively) 
[32] , and the fixed yields by ±50% of the theoretical val- 
ues [34] . To make a theory-independent determination of 
the impact due to T(lDj), the fit is also repeated with 
four of the T{lDj) 7X6j(lP) yields free to vary (the 
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FIG. 5: Fit to the 180 < E* < 300 MeV region of the T{3S) data (a) for all of the data, and (b) after subtraction of the fitted 
background contribution. The jndof (where ndof stands for number of degress of freedom) for the fit is 119.1/110. The thin 
gray lines indicate the individual signal components in the fit, as labelled. 
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TABLE II: Summary of the analysis of the 180 < E* < 300 MeV region of the Y{3S) data. The E* column lists the transition 
energy assumed in this analysis. Errors on the yield are statistical only. Regarding the derived B(x6j(2P) — > ^T{2S)): the 
BABAR value is from this paper, while the CUSB and CLEO columns are derivations based on [10] and [12] using up-to-date 
secondary branching fractions from [27]. For the BABAR result, the listed uncertainties are statistical, systematic, and from the 
uncertainties on secondary branching fractions, respectively. For the other results, the total uncertainty (all sources combined 
in quadrature) is given. Upper limits are given at the 90% confidence level. 

Transition E* Yield e Derived Branching Fraction (%) 

(MeV) (%) BABAR CUSB CLEO 

X6o(2P) ^ 77(25) 205.0 -347 ± 209 0.105 -4.7 ± 2.8lLo l ± 0-5 (< 2.8) 3.6 ± 1.6 < 5.2 
Xbi {2P) ^ -fT{2S) 229.7 4294 ±251 0.152 18.9 ± 1.1 ± 1.2 ± 1.8 13.6 ± 2.4 21.1 ± 4.5 
Xb2{2P) ^ -/T{2S) 242.3 2462 ± 243 0.190 8.3 ± 0.8 ± 0.6 ± 1.0 10.9 ±2.2 9.9 ±2.7 



T{lDi) ^ 7X6i(l^) and TilD^) ^ 7Xb2(lP) yields arc 
fit as a single component because their E* values are 
nearly identical, and the T(lDi) — > 7X&2(1^) transition 
is overvirhelmed by the main Xbi,2{'^P) jT{2S) peaks 
and remains fixed). Under this scenario, none of the 
r(lZ)j)-related transitions is found to be significant, and 
the yields are consistent with the theoretical predictions 
within statistical uncertainty. The Xb.jC^P) jT{2S) 
yields are not significantly affected. The changes in the 
fit yields for all of these alternative cases are added in 
quadrature and taken as the systematic uncertainty due 
to T{lDj) decays (^ 2%). It is worth reiterating that 
the excellent resolution obtained by using converted pho- 
tons separates the Y{lDj)- and X(,j(2P)-related compo- 
nents in E*, which is why the impact of the T(lDj) 
states does not dominate the measurement uncertainty. 
The fit is repeated with the bottomonium masses (hence, 
E* values) varied according to the PDG uncertainties 
[27], and the change in the yield added in quadrature 
(~ 2%). The number of Y(mS) mesons and its uncer- 
tainty (1.0%) were calculated separately, based on visible 
cross sections computed from dedicated e+e~ e+e~ (7) 
and e~^e~' — > /i+/i^(7) control samples. Systematic ef- 
fects due to the fit mechanics were tested by repeating 
the fit separately with an expanded E* range and a bin 
width of 0.5 MeV, the difference in results defining a small 
systematic uncertainty (1.5%). As a cross-check, the fit 
was repeated with the xtoC^P) component restricted to 
a physical range. The effect on the other signal yields 
was found to be small (< 2%). Finally, the background 
shape was replaced by a fifth-order polynomial and half 
of the resulting change in the yield (< 1%) taken as the 
symmetric error due to this assumed parameterization. 

We find B(r{'3S) ^ 7X5,7(2P)) x B{xbj ^ 7^(25)) = 
(-0.3±0.2l^ j)%, (2.4±0.1±0.2)%, and (1.1±0.1±0.1)% 
for J = 0, 1, and 2, respectively. Using B{T{3S) 
-/Xbj{2P)) from [27], we derive S(xw(2P) ^ 7^(2^)) = 
(-4.7±2.8l[]:^±0.5)%, (18.9±l.l±1.2±1.8)%,and (8.3± 
0.8±0.6± 1.0)%, where the errors are statistical, system- 
atic, and from the imcertainty on B{T{3S) — > 7X5,7(2P)), 
respectively. From these values, we calculate a 90% confi- 
dence level upper limit of S(x6o(2P) -)■ 7^(25)) < 2.8% 
[35]. Past experimental results [10, 12] averaged by the 
PDG [27] rely on assumptions for the branching fractions 



of r{2S) (+(r and T{2,S) 7Xb,/(2P) and their un- 
certainties that are no longer valid. In Table II, we have 
rescaled these previous results using the current values in 
order to make a useful comparison. We find our results 
to be in good agreement with the previous results, and 
to be the most precise values to date for the J =1 and 2 
decays. 



V. Y{fiS) : 300 < < 600 MeV 

The 300 <E*< 600 MeV range in the inclusive T(3S') 
photon energy spectrum, shown in Fig. 6, is compli- 
cated by many radiative bottomonium transitions. A 
principal feature is the photon lines from the three di- 
rect T{2>S) — >■ 7Xbj(lP) decays. Photons from the sec- 
ondary decays, X6j(l-P) — ^ I'^i^S), have energies that 
overlap with these initial transitions. There are sev- 
eral ways to produce Xbj(l-P) from T{2,S), each with 
unique Doppler broadening and relative rate. These de- 
cays "feed-down" to produce many extraneous Xb.ii^P) 
mesons that contribute substantially to the background 
level through subsequent Xbj{^P) I'^i^S) decay. At 
the lower edge of this energy range, there are potential 
contributions from T(3S') — > -fr]b{2S) and T{2S) produc- 
tion from initial state radiation (ISR). 

The best known of the T{3S) — >■ 7X67(1P) branching 
fractions comes from the CLEO experiment, which was 
able to isolate the T{3S) 7x60 signal [14]. A sep- 
arate analysis of X6j(lP) decays to multihadronic final 
states further set upper limits on B{T{3S) 7X6j(lP)) 
[36]. A recent analysis of T{3S) -> 7X61.2 (1^") 
77T(1S') transitions with exclusive Y(IS) — >■ £'^i~ 
decays has resulted in a measurement of T(3S) 
7X61, 2(1^*) branching fractions [15]. Our improved E* 
resolution with the converted photon sample allows us to 
disentangle the overlapping photon lines to make a direct 
measurement of these radiative transitions as well. We 
also search for a signal for T{3S) — >■ jr]b{2S). 

The direct T{3S) — > 7X6j(lP) lineshapes are parame- 
terized using the double-sided Crystal Ball function de- 
scribed in Sec. IV plus an independent Gaussian to ac- 
count for broadening from non-linearities in the E* reso- 
lution due to low momentum tracks encountered in this 
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energy range. The T(35) — > 7775(25') lineshape is mod- 
eled with the convolution of a rclativistic Breit-Wigner 
function (natural lineshape for the rii,{2S)) and a Crystal 
Ball function (experimental resolution function), where 
the Breit-Wigner function has been modified by a trans- 
formation of variables to E* using Eq. (1). The ISR- 
produced T{2S) signal is parameterized with a Crystal 
Ball function, for which the width is dominated by the 
spread in the e+e~ beam energy. 

The lineshapes for the decays Xbj{^P) — >■ 7^(1<S') de- 
pend on the initial decays that produced the Xbj{^P) 
states. We consider six main production pathways: 

. r(35) ^ 7X67(1P) 

• r(35) ^ 7ir{2S) ^ 777X6j(lP) 

. r{3S) ^ jjr{lDj) ^ 777X67(1P) 

• r(35) 7r7rr(25) nn-fXbji'i-P) 

• r{3S) ^ 7X67(2P) ^ jnnXbjilP) 

• e+e" 'yisRT{2S) lisRlXbj{^P)- 

The feed-down contribution from T(35) — >■ 7Xf,j(lP) 
is determined directly from the fit to the data. The line- 
shapes for the subsequent Xbj{^P) l^i^S) decays are 
distorted by Doppler-broadening effects. Wo parameter- 
ize the Xf,j(lP) transition lineshape with the convolution 
of a rectangular function and a Crystal Ball function. 
Because of the large Doppler width (~ 20MeV), the re- 
sulting shape is relatively broad and non-peaking. In the 
fit, the relative yields of the direct to the secondary tran- 
sitions are fixed according to the ratios of the expected 
efficiencies for each mode, and the branching fractions for 
the Xbj{^P) — >■ 7^(1'S') decays (to be discussed below). 

There are two 37 pathways from T(35) to x;,j(lP). 
Decays via T{2S) are fairly well understood, and the 
precision branching fraction results from Sec. IV are 
used to determine the expected yields and uncertainties. 
In contrast, the decays via T(lDj) have not been mea- 
sured in detail. We rely on theoretical predictions [34], 
found to be consistent with an experimental measure- 
ment of the 47 cascade to T{1S) [31], to estimate the to- 
tal feed-down component. We take the uncertainties on 
S(T(3S') 7X6j(2P)) [27] and introduce a 30% uncer- 
tainty on each theoretically calculated branching fraction 
in the decay chain. Doppler effects introduce a smooth 
~ 5 MeV broadening in these (and other) multi-step de- 
cay processes, thus the lineshapes for the individual 87 
pathways are adequately parameterized using a standard 
Crystal Ball function. 

There are two di-pion decay chains leading to Xbj(l-P): 
either via T{2,S) 'k-kT{2S) or X6j(2P) 7r7rXbj(lP). 
The former has been precisely measured by BABMt in 
a recent analysis of the recoil against tt'^tt^ to search 
for the /if, (IP) state [37]. We combine the branching 
fraction from that analysis with the PDG average [27] 
to obtain B{T{iS) iT+TT-r{2S) = (2.7 ± 0.2)%. For 
the 7r°7r° transition, we use the current world average 



branching fraction value [27]. The relevant MC sam- 
ples are generated with the experimentally-determined 
m^+TT- distribution [38]. Di-pion transitions between 
Xbj(2P) and Xbj(l-P) for J =1 and 2 have been measured 
experimentally by CLEO [39]. The above-mentioned 
BaBAR di-pion analysis [37] also measured these quanti- 
ties, which are averaged with the CLEO results to derive 
B(xbj(2P) 7r+7r-Xb.7(lP)) equal to (9.1 ± 1.0) x IQ-^ 
and (5.0 ±0.6) x 10^'^ for J 1 and 2, respectively. De- 
cays to the J = state, with different initial and final J 
values, and via tt^tt^ have thus far been below the level of 
experimental sensitivity. To calculate the expected feed- 
down, we assume isospin conservation such that T^a^a = 
\T^+^-, and estimate B{xbo{2P) ~^ -KTrxboi^P)) to be 
about one-fifth of that of the other J states [40]. We 
assume a 30% uncertainty on all theoretically-estimated 
branching fractions. 

Radiative decay of ISR-produced T{2S) mesons can 
yield Xbj(l-P) signals. The estimated production cross 
section for T{2S) is (28.6 ± 1.4) pb [41], where we have 
assigned a 5% uncertainty to this theoretical calculation. 
We combine this with the T{2S) 7Xfcj(l-P) branching 
fraction [27] to determine the size of this contribution to 
the background. From MC simulation, we conclude that 
the lineshape may be parameterized with a Crystal Ball 
function. 

Except for feed-down from T{3S) 7X6j(lP), which 
is determined from the data, the yields of these compo- 
nents are fixed in the fit. The branching fractions for 
the final step of the decay chain, ;B(x6j(lP) jT{lSj), 
are measured precisely for J = 1 and 2 in Sec. VI. Our 
values for these decays are averaged with results from 
CLEO [15]. For decays with J = 0, the CLEO [15] Col- 
laboration has recently presented observations. Since we 
do not observe this decay in Sec. VI, we use the measured 
branching fraction value from CLEO [15]. 

In the fit, we include two components related to 
hb{lP) — >■ jr]b{lS) decays. The hb{lP) decay is as- 
sumed to decay with a large branching fraction via 
/ib(lP) 7??h(lS') [42]. The two relevant hb{lP) pro- 
duction mechanisms are T{3S) — >■ 'K~^n~ hb{lP) and 
r(35) -)■ TT^hbilP). BaBar has studied both of these 
modes, finding B{r{3S) 7r+7r-/if,(lP)) < 2.5 x lO""^ 
[37] and BiT{3S) ^ n°hb{lP)) x B{hb{lP) lVb{lS)) 
= (4.7 ± 1.5 ± 0.6) X 10-4 [43]. Due to the effects 
of Doppler broadening, we parameterize the decay via 
7r° using the Doppler-broadened Crystal Ball function 
as described for Xbj{^P) 7^(15') transitions from 
T{3S) — >■ 7Xb7(l-P), and via tt+tt" using a standard 
Crystal Ball function. The yields for these components 
are fixed in the fit, and are nearly negligible. 

In the fit, all of the lineshape parameters are fixed 
to the MC-determined values except for the yield of the 
T{3S) lXbj{lP) (and its related Xbj{lP) ^ 7^(15") 
components), an overall E* scale offset, and the back- 
ground lineshape parameters. The feed-down yields are 
fixed using the branching fractions as described above. 
Repeated trials of the signal extraction on simulated 
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TABLE III: Summary of the analysis of the 300 < -B; < 600 MeV region of the r{3S) data. The E* column lists the transition 
energy assumed in this analysis. Errors on the yield arc statistical only. For the Derived Branching Fraction, the BABAR values 
are from this work, and the CLEO results are from [14, 15]. The upper limit is given at the 90% confidence level. 



Transition 


e; 

(MeV) 


Yield 


(%) 


Derived Branching Fraction (xlO ^) 
BABAR CLEO 


r{3S) - 
r{3S) - 

T{3S) - 


7Xi>2(lP) 
■>7Xw(lP) 
■> 7X60 (IP) 


433.1 
452.2 
483.5 


9699 ± 318 0.794 
483 ±315 0.818 
2273 ± 307 0.730 


10.5±0.3log 7.7 ±1.3 
0.5±0.3j:;5 ? (< 1.0) 1.6 ±0.5 
2.7 ±0.4 ±0.2 3.0 ±1.1 



datasets determine that, given the low efficiency and ex- 
pected number of events, and high level of background, 
obtaining a reliable yield for rn,{2S) and ISR-produced 
T{2S) is not possible. These components are there- 
fore not included in the fit. The measured photon en- 
ergy spectrum and the fitted yields are presented in 
Fig. 6, before and after the subtraction of the inclu- 
sive background. There is a clear separation of the 
T{?>S) — > "iXhj{^P) transitions, enabling us to observe 
the transitions to J = 0, 2, and find only a very small in- 
dication for J = 1. Table 111 summarizes the fit results. 

We consider systematic uncertainties due to the choice 
of background shape (1 — 2%), fit range, and bin- 
ning (1.5%), the effect of fixing parameters to the MC- 
determined values (3 — 6%), uncertainty in the photon 
conversion efficiency (3.6%), uncertainty in the T{rnS) 
counting (1.0%), uncertainty in the bottomonium masses 
(1 — 3%), and the impact of fixed feed-down yields 
(2%). The values in parentheses are representative of the 
T(3S') 7X60,2(1^') decays; for the Xbi-^lated results, 
the effects of the feed-down lincshapcs and the yields and 
the background shape dominate (about 20% each) due 
the marginal signal size. The evaluation of these uncer- 
tainties is done as described in Sec. IV, with the excep- 
tion of the feed-down-related uncertainty that is unique 
to this energy region. To assess the uncertainty related 
to the assumed branching fractions, we repeat the anal- 
ysis many times with the value of each input branching 
fraction varied randomly within its total uncertainty. We 
adopt the standard deviation of the change in the results 
as a systematic error. As a cross-check, we repeat the fit 
with the yields of the T(lZ?)-related feed-down compo- 
nents allowed to vary as a free parameter. We find only 
a small change (< 2%) in the overall branching fraction 
results, and consider this to be sufficiently accounted for 
by the systematic uncertainty determined from our pro- 
cedure of varying the branching fractions. Including ISR 
and r]b{2S) components in the fit produces an effect of 
less than ~ 1%, due to their slight impact on determin- 
ing the overall background shape. 

We measure B{T{3S) 7Xfej(lP)) = (2.7±0.4±0.2) x 
10-3, (0.5±0.3t!^;?) X 10-3, and (10.5 ± CStg;^) x 10-^ 
for J = 0, 1 and 2, respectively. We observe evidence for 
the T{3S) — >■ 7X60, 2(1^*) transitions, with total signifi- 
cances greater than 6.8a and 16a, respectively. We do 
not find evidence for the suppressed T{3S) 7Xfci(l-P) 
decay, and set the 90% confidence level upper limit of 



B(r(35) 7X61 (IP)) < 1.0 X 10-3. These results are 

consistent with previous limits [36], and improve upon 
the only measured value for the J = transition [14]. 
Our measurements of the T{3S) 7X61,2 (1-P) branch- 
ing fractions both differ from the recent CLEO observa- 
tions [15] by nearly 2a. Forcing the X6i,2(l-P) yields in 
our fit to match the CLEO results gives a poor /ndof 
of 399/293. However, using the 6(x6i,2(lP) ^ 7^(15")) 
results from Sec. VI to derive a total T(3S') "f"fT{\S) 
branching fraction via X6i.2(lP) (comparable to "J =1 
and 2" [15]), we find the results of the two experiments 
to be in close agreement. 

Adopting these results, we search for the T(3S') 
jr]b{2S) transition in the range 335 < E* < 375 MeV 
and find no evidence. Taking into account the domi- 
nant statistical uncertainty, we derive an upper limit of 
B{r{3S) 7?76(2S')) < 1.9 X 10-3 at the 90% confidence 
level. This limit is a factor of two larger than the limit 
set by CLEO [14]. 



VI. r{2S) : 300 <E*< 800 MeV 

We study five possible signals in the 300 < E* < 
800 MeV range in r{2S) data: three X6j(lP) ^ 
7T(1S') transitions, ISR r{lS) production, and T{2S) 
777fc(lS'). This energy region, shown in Fig. 7, has been 
analysed using calorimeter-detected photons by both 
BaBAR [18] and CLEO [19], the former finding evidence 
to confirm the r]b{lS). The improvement in resolution 
from the converted photon sample could allow a precise 
measurement of the r]b{lS) mass. However, because E* 
for the r{2S) -)■ -yVbilS) transition is « 613 MeV (com- 
pared to w 920 MeV in the T(35') data), its measure- 
ment is more difficult due to a lower detection efficiency 
and larger inclusive photon background. Studying this 
energy range is nonetheless useful, since the branching 
fractions for X6j(lP) 7^(15) have had large uncer- 
tainties [4, 6, 7] until very recently [15], and the values 
are necessary inputs to the analysis described in Sec. V. 
The J = decay has also only been recently observed 
[15? ]. These external measurements were unavailable 
when this analysis was initiated. 

We parameterize the X6j(l-P) transition lineshape with 
a Doppler-broadened Crystal Ball function, as described 
in Sec. V. The ISR and T{2S) 7r?b(15) lineshapes 
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FIG. 6: Fit to the 300 < E* < 600 MeV region of the T{3S) data (a) for all of the data, and (b) after subtraction of the 

fitted background contribution. The jndo] for the fit is 316/290. The thick solid lines indicate the total fit, the thin 

solid lines indicate the T(3S) 7x6 j (If*) components, and the dashed lines indicate those from X6j(lP) — >• 7T(1S) and 
hb{lP) — >■ 7r?6(lS'), as labeled. 
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are modeled with a Crystal Ball function, and relativis- 
tic Breit-Wigner function convolved with a Crystal Ball 
function, respectively. The Iinc;sliape parameters are de- 
termined from MC samples. Several different natural 
widths are tested for the T]t,{lS), and because the Crys- 
tal Ball parametcir values (related to E* resolution) are 
found to be independent of the width, the values averaged 
over all samples are used. In the fit to the data, all of 
the parameters arc fixed to these MC-determined values, 
except for the yields for the X6j(l-P)) ISR, and r]b{lS) 
signals, the mass of the r]b{lS), the inclusive background 
shape parameters, and an overall E* scale offset. The 
width of rjbilS) is fixed to 10 MeV. 

Figure 7 shows the converted photon energy spectrum 
before and after the subtraction of the inclusive back- 
ground, with an inset focusing on the region of the ex- 
pected T{2S) jijbilS) transition. The E* resolu- 
tion provides clear separation of the x;,i^2(l^')-i'elated 
peaks, allowing for the first direct measurement of these 
transitions in an inclusive sample. The results of the 
fit are summarized in Table IV. We find no evidence 
for X6o(l^') ^ lT{lS) decay. The r{lS) yield from 
ISR production is consistent, within large uncertainties, 
with the result scaled from the previous BaBAR mea- 
surement [18]. As expected from signal extraction stud- 
ies on simulated datasets, the search for a signal in the 
T{2S) — >■ 7?7f,(lS') energy region does not find a rehable 
result. Estimating the statistical significance from the 
change in of the fit with and without this component 
results in the equivalent of a less than 2.5a effect. The 
E* scale offset in this energy range is — MeV. 

The systematic uncertainties on these measurements 
are related to the choice of background shape, the fit 
mechanics, the effect of fixing parameters to the MC- 
determined values, uncertainty in the photon conversion 
efficiency, uncertainty in the T{mS) counting, uncertain- 
ties in the bottomonium masses, and assumptions on the 
r]b{lS) width. The methodology for the evaluation of 
these uncertainties has been described for the most part 
in Sec. IV. The systematic uncertainty related to the 
r]b{lS) width is estimated by finding the maximal change 
in yield when the fit is repeated using a range of widths 
between 2.5 — 15 MeV, values consistent with a wide range 
of theoretical predictions. While varying the assumed 
rib{lS) width affects the event yield, it is found to have a 
negligible impact on the significance of the signal. For the 
Xbi,2i^P) — > 7T(1S') transitions, the largest sources of 
uncertainty are related to the fixed lineshape parameters 
(3 — 4%), uncertainty in the bottomonium masses (~ 4% 
for Xf)2(2P), and dominant for the E* scale uncertainty) 
and the conversion efficiency (5.2%). Each of the remain- 
ing sources contributes less than 2%. For the rib{lS) sig- 
nal, systematic uncertainties dominate the result. The 
largest effects are due to varying the background shape 
(^ 31%), the bottomonium masses (~ 25%), the MC- 
determined parameters (~ 22%), and the r]b{'i-S) width 
(~ 16%). 

We measure B{T{2S) ^ 7X&j(l^)) x B{xbj{lP) ^ 



7T(1S')) = (8.3±5.6t^;^)xl0-4, (24.1±0.6±1.5) x lO^^, 
and (13.9 ± 0.5tH) x IQ-^, for J = 0, 1 and 2, respec- 
tively. Using 6(r(25) -)■ 7Xbj(lP)) from the PDG [27], 
we derive B{xbj{lP) -> 7^(15')) = (2.2±1.5lJ:^±0.2)%, 
(34.9 ± 0.8 ± 2.2 ± 2.0)%, and (19.5 ± O.TtH ± 1.0)%, 
where the uncertainties are statistical, systematic, and 
from the uncertainty on B{T{2S) jXbj{^P)), respec- 
tively. We calculate a 90% confidence level upper limit 
of Bixboi'^P) lT{lS)) < 4.6%. As previously, we 
rescale the existing results [4, 6] using the most up-to- 
date secondary branching fraction values [27] to obtain 
the results quoted in Table IV. Our X5j(l^') transition 
results agree with the previous measurements, but rep- 
resent a two- to three-fold reduction in the total uncer- 
tainty. We find reasonable agreement with, and a compa- 
rable precision to, the recent measurements from CLEO 
[15]. When the yield-related systematic uncertainties on 
the measurement of the r]b{lS) candidate are taken into 
account (excluding those due to the 776(15*) width), the 
result is further reduced in significance to an equivalent of 
~ 1.7(7. We find no evidence for an r]b{lS) signal in this 
analysis of the T{2S) dataset, and set a corresponding 
limit of B{r{2S) 'ymi'i-S)) < 0.21%. 

VII. r(3S) : 600 <E*< 1100 MeV 

The analysis of the 600 < E* < 1100 MeV region 

for the T{3S), shown in Fig. 8, is very similar to 
that in Sec. VI of the 300 < E* < 800 MeV region 
for the T{2S). Again, we study potential signals from 
three Xbj('2P) — > 7T'(1S') transitions, T{1S) production 
from ISR, and T{3S) ■^rjbi^S). In this case, the 
calorimeter-based analysis of the same region produced 
the discovery of the ?7f,(lS') [17]. The higher E* value for 
T{3S) 777h(lS') offers the advantages of both an in- 
creased efficiency and lower background level compared 
to the analagous analysis in T{2S) data, and therefore 
a better sensitivity for the observation of r]b{lS). There 
is also the possibility of updating the measurements of 
Xbj{^P) — >■ ^T{1S) transitions, including confirmation 
of the decay of the J = state [10, 12]. 

We parameterize the signal lineshape in the same man- 
ner as described in Sec. VI, with Doppler-broadened 
Crystal Ball functions for the Xbj{^P) transitions, a 
Crystal Ball function for ISR production of the T{1S), 
and the relativistic Breit-Wigner Crystal Ball convolu- 
tion for the rib{lS) signal. As before, all of the lineshape 
parameters are fixed to their MC-determined values, with 
the yields for the Xfcj(2-P), ISR, and rib{lS) signals, the 
mass of the rib{lS), the inclusive background shape pa- 
rameters, and an overall E* scale offset free to vary in 
the fit. An rib{lS) width of 10 MeV is assumed. 

Figure 8 shows the converted photon energy spec- 
trum and fitted yields before and after the subtraction 
of the inclusive background, with an inset focusing on 
the E* region of the expected T{3S) — > jrjb{lS) transi- 
tion. The results are summarized in Table V. Although 



14 




8 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 



(GeV) 



> 



2000 



1500 



^ 1000 



500 







_L I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I M 

:(b) 



% 300 

200 



g ")(> 



-100 

X (IP) 



ISR 




11 (IS) 



0.5 0.55 0.6 0.65 0.7 

Ej(GeV) 



ISR n(is) 



,(ip) 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I- 



0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 

E* (GeV) 



FIG. 7: Fit to the 300 < < 800 MeV region of the T{2S) data (a) for all of the data, and (b) after subtraction of the 
fitted background contribution, where the inset focuses on the T(25) — > 7r?6(15) region of the fit. The thin lines indicate the 
individual fit components. For this fit, /ndof = 511.0/487. 
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TABLE IV: Summary of the analysis of the 300 < E* < 800 MeV region of the Y{2S) data. The E* column hsts the transition 
energy assumed in this analysis, or in the case of T{2S) — > 7776(1<S'), the most significant (~ 1.7cr) feature in the relevant 
E* region. Errors on the yield are statistical only. Regarding the derived branching fractions, the BABAR value is from this 

paper, while the Crystal Ball (CB) and CUSB columns are derivations based on Ref. [4] and [6] using up-to-date secondary 
branching fractions from the PDG [27]; the CLEO results are from [15]. For the BABAR result, the listed uncertainties are 
statistical, systematic, and from the uncertainties on secondary branching fractions, respectively. Upper limits are given at the 
90% confidence level. Dashes indicate that no value has been reported in the relevant reference. 

Transition E* Yield e Derived Branching Fraction (%) 

(MeV) (%) BABAR CB CUSB CLEO 

Xbo{iP) 7^(15') 391.5 391 ± 267 0.496 2.2 ± 1.5^^^ ± 0.2 (< 4.6) <5 <l2 1.7 ±0.4 
XM (IP) ^ 7T(1S') 423.0 12604 ± 285 0.548 34.9 ± 0.8 ± 2.2 ± 2.0 34 ± 7 40 ± 10 33.0 ± 2.6 
Xw(lP) ^ 7T(1S') 442.0 7665t 272 0.576 1 9.5 ± 0.7li-^ ± 1.0 25 ± 6 19 ± 8 18.5 ± 1.4 
r{2S) ^ ■yVijlS) eiS.ltl llli 1109 ±348 1.050 0.11 ± 0.04tg g^ (< 0.21) - 



the X6i,2(l-P)-related peaks overlap, the E* resolution 
is still sufficient to measure the separate contributions. 
We find no evidence for Xbo{'2P) ■jT{lS) decay. The 
T{1S) yield from ISR production is in greement with the 
expectation from the previous BABAR measurement [17]. 
The best fit for a signal in the E* range corresponding 
to T{3S) 7776(15) has E* « 908 McV, which is a de- 
parture from, but not significantly inconsistent with, the 
nominal PDG value of 920.6^|^MeV [27]. Estimating 
the statistical significance from the change in x^ of the 
fit with and without this component results in the equiv- 
alent of a less than 2.9a effect. Based predominantly on 
the positions of the Xf)i.2(2P) transition peaks, the E* 

scale offset in this energy range is — O.Olg g MeV. We fur- 
ther verify that the E* scale is correct by repeating the fit 
with the peak positions of the Xbj{2P) and ISR compo- 
nents allowed to vary, and they are found at the expected 
locations. We also repeat the analysis with the E* scale 
offset forced to reproduce an %(15') result corresponding 
to the E* value for the nominal m^^(i5). The assumption 
that the observed mass difference is due to an offset in 
the energy scale by ~ 12 MeV is completely inconsistent 
with the photon energies observed for the well-established 
X6i.2(2P) states. Even with only a 5 MeV shift, the fit 
returns X6j(2P) — >■ jT{lS) yields that disagree with the 
world average [27] by more than a factor of two, and a 
X^/ndof « 840/492. 

The sources of systematic uncertainty and their evalua- 
tion are identical to those listed in Sec. VI. The main dif- 
ference between the two energy regions is that, as previ- 
ously remarked, the improved efficiency and background 
conditions in the 600 < E* < 1100 MeV region of the 
T(3S') dataset lead to fit results that are more stable. 
For the X6i,2(2-P)-related measurements, the dominant 
systematic uncertainty is due to the conversion efficiency 
(3.6%), and aU other sources are less than 2%. For the 
rib{lS) signal, the largest uncertainty in the yield is re- 
lated to the assumed 775(15) width (^27*^)- Of the re- 
maining systematic uncertainties, the largest two are due 
to the MC parameterization (~ 15%) and bottomonium 
masses (~ 4%), both enhancing the yield in a positive 
direction. Uncertainty due to the background shape. 



the largest factor in the equivalent T(25) analysis, is 
well controlled in the T(35) dataset and contributes less 
than 3% to the total uncertainty. The uncertainty in E* 
is dominated by statistical uncertainty, and the largest 
systematic contribution is related to uncertainty in the 
E* scale via the uncertainty in the other bottomonium 
masses [27]. 

We measure S(T(35) ^ 7Xbj(2P)) x fi(xbj(2P) ^ 
7r(15)) = (3.9±2.2tjj) X 10-4, (12.4±0.3±0.6) x IQ-^, 
and (9.2 ± 0.3 ± 0.4) x lO'^, for J = 0, 1 and 2, respec- 
tively. Using S(r(35) 7Xbj(2P)) from the PDG [27], 
we derive S(xfcj(2P) ^ 7^(1^)) = (0.7±0.4t°:?±0.1)%, 
(9.9±0.3t!^j±0.9)%, and (7.0 ±0.2 ±0.3 ±0.9)%, where 
the uncertainties are statistical, systematic, and from 
the uncertainty on B(T(35) — >■ 7Xbj(2P)), respectively. 
From these values, we calculate a 90% confidence level 
upper limit of S(xfco(2P) ^ 7^(15)) < 1.2%. As before, 
we rescale the previous results [10, 12] using the relevant 
branching fractions [27] to produce the values for com- 
parison in Table V. For the S(x6o(2P) 7^(15)) value 
from CUSB II [10], we convert the result to an upper 
limit of < 1.9% at the 90% confidence level. Our Xbj(2P) 
transition results agree with the previous measurements, 
and arc the most precise measurements to date. Assum- 
ing the peak near E* = 900 MeV to be due to decays 
to r76(15), our best fit result is B(r(35) -)■ 77/6(15)) = 
(5.8 ± 1.6tl]t} X lO"'*- The total significance of this re- 
sult once systematic uncertainties are included is ^ 2.7fT, 
and we set a limit of B{T{3S) 7?76(15) < 8.5 x lO"". 
We repeat the fit with the T]b{lS) mass constrained to 
the PDG value and its uncertainty [27]. The signifi- 
cance of this constrained result is < 1.9a. Wc measure 
B{T{3S) 77?6(15) = 3.8 ± which translates 

into an upper limit of B(r(35) ^ 7^/6(15) < 6.1 x IQ-^. 

VIII. DISCUSSION 

To conclude, we review the results of this study and 
their broader implications. The results for B{xbj{'nP) — > 

7T(?77,5)) presented here are the first derived directly 
from a measurement of the photon spectrum. For J = 1 
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TABLE V: Summary of the analysis of the 600 < E* < 1100 MeV region of the r{3S) data. The E* column hsts the transition 
energy assumed in this analysis, or in the case of T{3S) — ^ 7776(16'), the most significant (~ 2.7a) feature in the relevant E* 
region. Errors on the yield are statistical only. Regarding the derived branching fractions: the BABAR value is from this paper, 
while the CUSB and CLEO columns are derivations based on [10] and [12] using up-to-date secondary branching fractions 
from [27]. For the BABAR result, the listed uncertainties are statistical, systematic, and from the uncertainties on secondary 
branching fractions, respectively. Upper limits are given at the 90% confidence level. Dashes indicate no value has been 
measured in the quoted reference. 



Transition 



e; 

(MeV) 



Yield e Derived Branching Fraction (%) 

(%) BABAR CUSB CLEO 



Xw(2P) ^7^(15) 
Xti{2P) ^ jT{lS) 
X62(2P)^ 7^(15) 

r(35) ^ -rvbjis) 



742.7 
764.1 
776.4 
907.9 ±2.8 ±0.9 



'-269 

149651^^3 1-039 
11283tit 

j+263 
-262 



4693 



933: 



1.025 0.7±0.4t^;^ ±0.1 (< 1.2) 



9.9 ± 0.3 



-1-0.5 
-0.4 



±0.9 



1.056 7.0 ±0.2 ±0.3 ±0.9 
1.388 0.058 ± 0.016t°o°o\e (< 0.085) 



< 1.9 < 2.2 
7.5 ±1.3 10.4 ±2.4 
6.1 ±1.2 7.7 ±2.0 



TABLE VI: Comparison of the experimental branching frac- 
tion results from this work (BABAR) and some theoretical pre- 
dictions [34]. 



Decay 


BABAR. (%) Theory (%) 


B(Xi,o(2P)^ 7^(25)) 


(< 2.8) 


1.27 


B(X6i(2P)^7T(25)) 


18.9 ±2.4 


20.2 


B(X62(2P)^ 7^(25)) 


8.3 ±1.4 


10.1 


Bixto(2P) ^ yniS)) 


« 1-2) 


0.96 


B{xbi{2P) ^ jTilS)) 


9.9 ± 1.1 


11.8 


B{xb2{2P) ^ yTilS)) 


7.0 ± 1.0 


5.3 


B(X6o(lP) ^ 7^(15)) 


(< 4.6) 


3.2 


B(X6i(lP)^ 7^(15)) 


34.9 ±3.1 


46.1 


B(X62(1P)^ 7^(15)) 


19.5111 


22.2 



TABLE VII: Comparison of our results with predictions [44, 
45] for T{3S) — > IXtji^P) decays. We convert our result into 
partial widths (in units of eV) using a total width of rr(3S) = 
20.32 ± 1.85 keV [27], absorbing this additional uncertainty 
into the total. 



Source 


J = 


J = 1 


J = 2 


BABAR 


55 ±10 


< 22 


216 ± 25 


Moxhay-Rosner 


25 


25 


150 


Crotch et al. 


114 


3.4 


194 


D aghighian- Silverman 


16 


100 


650 


Fulcher 


10 


20 


30 


Lahde 


150 


110 


40 


Ebert et al. 


27 


67 


97 



and 2, we have made some of the most precise measure- 
ments of these branching fractions to date, thus helping 
to resolve some discrepancies between previous experi- 
mental results {i.e. in xtjC^P) jT{2S) decays). Table 
VI shows a comparison of our results with some theoret- 
ical predictions [34] . These predictions are in reasonable 
agreement with our experimental results. 

Our observations of T{3S) — >■ 7x50,2 (1-P) decays con- 
firm the general features seen in previous measurements 
[14, 15, 36]: decays to J = 1 are suppressed compared 
to J = 2 and 0. This is unusual compared to all other 
S ^ P radiative transitions in the heavy quarkonium 
system measured thus far. As noted previously [44], the 
wavcfunction overlap in the (3^51 | r | l^Pj) matrix 
elements is unusually small. Therefore, predictions for 
these decay rates are largely dependent on higher-order 
relativistic corrections and are thus sensitive to specific 
details of the chosen theoretical model. That said, the 
comparison of our results with a selection of theoretical 
predictions [44, 45] shown in Table VII (where we have 
converted our branching fraction measurements into par- 
tial widths) finds no good agreement with any particular 
model. Indeed, even the hierarchy of the decay rates 
( J = 2 > > 1) is generally not well predicted. Further 
work, both theoretical and experimental, will be required 
to understand these decays. 



The searches for r]b{lS) and r]b{2S) states using the 
converted photon energy spectrum are largely inconclu- 
sive. Over a range of approximately 9974 < m^^(2S) < 
10015 MeV/c^, we find B(r(35) -^77?6(25)) < 1.9x10-3. 
This value is consistent with, but does not improve upon, 
previous measurements [14]. Due to low efficiency and 
high background, no evidence for Y{2S) — >■ 'yrii,{lS) is 
found. In the T{3S) system, the most significant peak- 
ing structure in the E* energy region expected for the 
r{3S) -)• 77/6(15') transition has E* « 908 MeV. If inter- 
preted as an rjbilS) signal, this value trends toward the 
most recent potential model [46] and lattice [47] predic- 
tions, but we caution that the significance of this result 
is insufficient to draw such a conclusion regarding the 
r]b{lS) mass. Taking advantage of the improved resolu- 
tion from a converted photon technique to make a defini- 
tive measurement of the r]b{lS) mass and width will re- 
quire more data from future experiments. 
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TABLE VIII: Selection criteria for the e^e 7 and n 7 
efficiency studies. 

Quantity e^e~(p() /i"'"/i~(7) 
nTRK ^4 ^4 
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X2 > 0.50 > 0.75 

Greater I cose* I < 0.70 < 0.70 

Lessor I cos (9* I < 0.65 < 0.65 

o' (") < < 20 
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Appendix A: SYSTEMATIC UNCERTAINTIES 
ON MC-DETERMINED EFFICIENCIES 

Branching fraction measurements in this analysis rely 
on MC-generated signal decays to determine the photon 

conversion and reconstruction efficiency. This efficiency 
is dependent on the detector material model. To evaluate 
a systematic effect due to the understanding of the detec- 
tor in the simulation, a comparison of e+e~7 and /i+/i~7 
samples between data and MC is made. Inclusive de- 
cays to an e+e~ or /x^/z" pair plus a photon are selected 
by requiring exactly four charged tracks in the event. 
The CM momentum of the two highest-momentum non- 
conversion tracks as a fraction of half of the CM beam 
energy (xl, Xg), the higher and lower values of their 
CM polar angles ( | cos f?* 2 1 ) , and the CM acolinearity 
(a*), are used as discriminating variables. We require 
e"'"e~7 events to pass a predefined filter optimized to se- 
lect Bhabha scattering events, and for the /z+/Li~7 events 
to fail this requirement. In cases of multiple candidates 
per event, the candidate with m^+£-^ closest to the CM 
beam energy is retained. The values for the selection 



criteria variables are summarized in Table Vlll. 

To avoid contamination from resonant decays {e.g. 
Xbj{nP) 'yr{mS){e+i-), or r{nS) plus 
an extraneous photon), only the off-peak datasets are 
used for this study. The e~^e~{'-f) MC sample uses the 
BHWIDE generator [48], while the iX^n~^ MC sam- 
ple is generated using the KK2f generator [49]. The 
acceptance-based cross sections for these processes used 
in the MC generation arc calculated separately from this 
analysis as part of standard luminosity measurements in 
BABm. 

A systematic correction to the MC-dctermined effi- 
ciency is determined by comparing the number of events 
expected from the luminosity- weighted MC samples with 
the total number reconstructed in the data. The uncer- 
tainty on this correction (dominantly statistical) is used 
as the systematic uncertainty in the efficiency due to the 
detector material model. The four samples (e+e^7 and 
/i+/i^7 in off-peak T{2S) and T(3S') data) are averaged 
to calculate this number, as is justified by verifying excel- 
lent data-to-MC agreement across all relevant cos^, E*, 
and ranges. Integrated over all events, the ratio of the 
data and MC is 96.3 ± 3.1% when modeling the photons 
converted in the detector material. This value is applied 
as a correction factor, with 3.3% (when considering cross 
section uncertainties of about 0.8%) taken as an estimate 
for the systematic uncertainty in the efficiency. 

The MC-based signal efficiencies are also dependent 
on assumptions regarding inclusive bottomonium decays. 
The nTRK requirements attempt to select multihadronic 
final states. A difference in nTRK distributions between 
simulation and data could lead to an error on the recon- 
struction efSciency. To determine the size of this effect, 
the analysis is repeated with the requirements nTRK 
greater than 5 or nTRK greater than 6. The largest 
change in the efficiency-corrected yields for the most sig- 
nificant transitions (x6i,2(l,2P) ')T{1,2S)) is found 
to be 1.0%. 

Uncertainty in the modelling of the 7r° veto efficiency 

is tested in a similar manner, by repeating the analy- 
sis with the veto excluded and examining the change in 
the fit results for the most statistically significant tran- 
sitions. We rescale the MC-derived efficiency to equal 
half of the difference between the weighted average of 
the nominal and non-7r° vetoed results, and introduce a 
systematic uncertainty large enough to cover this differ- 
ence. Because lower energy photons are more suscep- 
tible to the application of a tt*^ veto, we find the dif- 
ferences to be energy-dependent and assign a different 
correction and uncertainty for each energy region. The 
scale factors (uncertainties) range from 0.991 (0.9%) for 
the 600 <E* < 1100 MeV range in T(3S') data to 0.963 
(3.9%) in the r(25) data. 

We combine these values to estimate a total systematic 
uncertainty on the efficiency of 3.6% to 5.2%, depending 
on the transition. 
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